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ABSTRACT 

The polytropic equation of state of an atomic hydrogen gas is examined for primordial halos 
with baryonic masses of Mh ~ 10 — 10 9 M©. For roughly isothermal collapse around 10 4 K, 
we find that line trapping of Lyman a (HI and Hell) photons causes the polytropic exponent 
to stiffen to values significantly above unity. Under the assumptions of zero H2 abundance and 
very modest pollution by metals (< 10~ 4 Solar), fragmentation is likely to be inhibited for such 
an equation of state. We argue on purely thermodynamic grounds that a single black hole of 
~ 0.02 — 0.003Afh can form at the center of a halo for z — 10 — 20 when the free-fall time is 
less than the time needed for a resonantly scattered Lyman a photon to escape from the halo. 
The absence of H 2 follows naturally from the high, > 10 4 K, temperatures that are attained 
when Lyman a photons are trapped in the dense and massive halos that we consider. An H 2 
dissociating UV background is needed if positive feedback effects on H 2 formation from X-rays 
occur. The black hole to baryon mass fraction is suggestively close to what is required for these 
intermediate mass black holes, of mass Mbh ~ 10 4 — 10 6 M , to act as seeds for forming the 
supermassive black holes of mass ~ 0.001M sp h e roid found in galaxies today. 

Subject headings: cosmology: theory - galaxies: black holes -ISM: clouds - ISM: atoms - atomic processes 
- radiative transfer 



1. Introduction 

A fundamental issue in the study of galaxy evo- 
lution is the formation of the central (supermas- 
sive) black hole. Accretion onto these black holes 
provides the energy source for active galactic nu- 
clei, which in turn impact the evolution of galaxies 
(Silk 2005). Earlier attempts at providing seeds 
for galactic black holes include dynamical friction 
and collision processes in dense young stellar clus- 
ters (Portegies Zwart et al. 2004) and formation 
from low angular momentum material in primor- 
dial disks (Koushiappas, Bullock & Dekel 2004). 

In this work, we consider the impact of the 
polytropic equation of state (EOS) of a metal- 
free, atomic hydrogen gas on the expected collapse 
of matter inside massive halos. The impact of a 
Solar metallicity polytropic EOS on the expected 



masses of stars in local galaxies has been inves- 
tigated by Li, Mac Low & Klessen (2003). The 
influence of molecular hydrogen and metal-poor 
environments has received detailed interest from, 
e.g., Abel, Bryan & Norman (2002, 2000); Bromm, 
Coppi & Larson (2002) for the formation of the 
first stars and from Scalo & Biswas (2002) and 
Spaans & Silk (2005) for the properties of the poly- 
tropic equation of state. From the work of Li et al. 
(2003) it has become clear that a polytropic EOS 
of state, P (x p 1 with p the mass density and 7 the 
polytropic exponent, strongly suppresses fragmen- 
tation of interstellar gas clouds if 7 > 1. This pa- 
per concentrates on the impact of Lyman a photon 
trapping on the EOS and the interested reader is 
referred to Rees & Ostriker (1977) and Silk (1977) 
for some of the fundamental thermodynamic and 
star formation considerations that come into play 



1 



here. 

2. Model Description 

We assume a metal-free hydrogen gas that is 
cooled by Lyman a emission as it collapses inside 
a dark matter halo and radiates away about twice 
its binding energy (Haiman, Spaans & Quataert 
2000). Note that Lyman a cooling is expected to 
dominate over radial contraction factors of at least 
15-60 as long as the metallicity is less than 0.1 of 
Solar (Haiman et al. 2000) . The absence of any H2 , 
that would cool the gas to below 8,000 K, is crucial 
in this and we come back to this point in Section 4. 
We further employ a polytropic EOS and a perfect 
gas law, P oc pT, for the gas temperature T, and 
write 7 as 



This last step is justified (Scalo & Biswas 2002) as 
long as the heating and cooling terms in the fluid 
energy equation can adjust to balance each other 
on a time-scale shorter than the time-scale of the 
gas dynamics (i.e., local thermal equilibrium). Be- 
low, we compute the polytropic EOS for the case 
that the cooling time is shorter than the free-fall 
time and for the case that the photon propagation 
time exceeds the dynamical time. 

It should be noted that because 7 depends on 
the (logarithmic) derivative of the temperature 
with respect to density, it implicitly depends on 
radiative transfer effects and changes in chemi- 
cal composition through derivatives of the heat- 
ing and cooling functions (Spaans & Silk 2000). 
The Lyman a radiative transfer techniques as de- 
scribed in Haiman & Spaans (1999) and Dijkstra 
et al. (2006) are used to compute the transfer of 
Lyman a photons. 

We consider spherical dark matter halos that 
have decoupled from the Hubble flow and are char- 
acterized by a mean density of p « 200/?b(l + z) 3 
at z = 10 — 20, for a baryonic number density 
Pb/iriH = 3 x 10~ 7 cm~ 3 today, hydrogen mass 
ma, total halo masses of Mf t = 10 7 — 10 9 M and 
a characteristic size scale of L = (SMh/Airp) 1 / 3 . 
This yields, over z = 10 — 20, a typical mean 
density and column of n = 0.05[(1 + z)/10] 3 
cm- 3 and N Q = 10 22 [(1 + z)/10] 2 (M, l /10 9 M Q ) 1 /3 
cm~ 2 , respectively. We further assume that mat- 
ter inside the halo remains at approximately 10 4 



K during its collapse so that an isothermal den- 
sity profile, n oc no(L/r) 2 , is applicable for every 
radius r. Therefore, the column a Lyman a pho- 
ton has to traverse from a radius r to L scales 
as J ndr ~ L/r — 1 with a mass inside of r of 
M(r) ~ r. 

In the absence of any ionizing sources, heat- 
ing is provided by gravitational compression, T oc 
n 1 - 5 . The velocity dispersion of the gas is ther- 
mal and equal to Ay = I2.9T4 km/s, with T4 in 
units of 10 4 K. The natural to thermal line width 
of the Lyman a line is denoted by a and equal to 
a = 4.7x 10- 4 T 4 ~ 1/2 . 

3. Results 

3.1. Static Case 

With cooling provided by Lyman a emission 
only, the thermal equilibrium of the baryonic 
matter in the halo approximately (Spitzer 1978; 
Haiman et al. 2000) follows, for r' in units of L, 

7.3 x 10- 19 n e (r')n H (r')e- 11S ' im / T ^e(r') = 

1.9n(r')GM h /Ltjj, (2a) 

with M h /L = M h (r'/L)/r' = M(r')/r', electron 
density n e , atomic hydrogen density uh, Lyman 
a escape fraction e and free-fall time 

t ff = 4.3 x 10 7 /™(r') 1/2 yr- (26) 

It is assumed here that the cooling time at the 
peak of the cooling curve is t c = 3/2nkT V i r /n 2 A, 
with Boltzmann's constant k, the halo's virial tem- 
perature T V i r and A w 2 x 10~ 22 erg s _1 cm 3 . De- 
pending on the ambient conditions, the medium 
cools somewhere around the peak of the cooling 
curve and T mr < 10 4 K if the mass is smaller than 
10 8 [(1 + z)/10]- 1 - 5 Mq for the virialization red- 
shift z (Haiman, Rees & Loeb 1997). 

The escape fraction e of Lyman a photons from 
a sphere diminishes from unity when collisional 
de-excitation above a critical HI column density 
iV c becomes important (Neufeld 1990). This col- 
umn N c depends on the ambient temperature and 
ionization balance through the probability for col- 
lisional de-excitation po = gp "^^ e " e , with pro- 
ton density n p , collisional de-excitation rate coef- 
ficients q p and q e , and the Einstein A coefficient 
A 2 i connecting the 2p and 2s states. In this, the 
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ambient proton density is assumed to be lower 
than about 10 4 cm~ 3 so that the created 2s hy- 
drogen atoms undergo two-quantum decay to the 
ground state. It follows that iV c ranges between 
10 21 and 10 23 cm" 2 for y = ripT^ - 17 between 10 2 
and 10 4 cm -3 , respectively (Dijkstra et al. 2006), 
and is much larger for much smaller proton densi- 
ties. These values for N c are a factor of a few larger 
than the corresponding values for a slab (Neufeld 
1990), since resonantly scattered photons escape 
more easily from a sphere than from a slab for the 
same surface-to-center optical depth. 

Furthermore, following the Monte Carlo radia- 
tive transfer techniques in Dijkstra et al. (2006) 
and Haiman & Spaans (1999), and for a HI col- 
umn N H , we find that e w N c N^-°, for N H = 
2N C — 100_/V C and for spherical clouds. In deriving 
this fit to the numerical results, we have made sure 
that the line profile is sampled far enough into the 
wings to accurately determine e and N c . When 
applied to a slab, rather than a sphere (see the 
analytical solution in the appendix of Dijkstra et 
al. 2006), our method yields results that agree well 
with those of Neufeld (1990, his Figure 18). 

For the resonantly scattered Lyman a line, it 
follows, for a line center optical depth To, mean 
line opacity a s and profile function 4>{x) in nor- 
malized frequency units x, that an escaping pho- 
ton, which scatters N times, experiences a fre- 
quency shift x s ~ iV 1 / 2 and travels a distance 
N 1 / 2 /(a s (j)(x s )) that is equal to the size of the 
medium To/a 8 . Hence, x s ~ to</>(x s ) ~ (aro) 1 / 3 
since <f> ~ a/irx 2 . On average, a time 5t ~ 
^/(to4>(x s )) elapses between the ~ N scatterings. 
Thus, a time t p h ~ NSt ~ -^(aTo) 1 / 3 is required 
for a photon to escape, where the optical depth 
is given by r = 1.04 x 10~ 13 N H T^ 1/2 . Typi- 
cally, we have that To > 10 7 . Thus, for a given 
density, and in the limit that tff >> t p h with 
Nh > N c , an increase in column leads to a pro- 
portional decrease in the spherical escape proba- 
bility 1 — e~ aT °/aT « 1/aro and the chance that a 
scattering hydrogen atom will not suffer collisional 
de-excitation effects. 

Obviously then, Equation (1) has a weak de- 
pendence of 7 on density for modest columns due 
to the exponential temperature dependence and 
the Nfl 1 ' cx n 2 ^ 3 scaling of e in the static case 
and for fixed mass Mh- Under collisional ioniza- 



tion equilibrium, one finds from solving Equations 
(2) for T(n) that 
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2log(Cn 1 / 2 ) 



= 0.006 - 0.007, 



(3) 



for proton densities larger than 10 cm , Nh > 
10 2i-23 cm -2 or r < rstat = i q _ Q.oiL for all 

halos over z = 10 - 20, and where C <~ 10~ 36 
Mh/10 7 M & cm 3 / 2 . Hence, as to be expected, the 
stiffening of the polytropic EOS is always modest 
when the exponential temperature dependence of 
the Lyman a cooling rate acts unchecked. 

3.2. Dynamic Case: HI 

For the halos considered here, t c < tff by a 
factor of a few. However, if the random walk that a 
Lyman a photon performs takes a time t p h that is 
comparable to or longer than the dynamical time 
on which the halo evolves, cooling is effectively 
shut down. Photons can then only escape through 
parts of the line wings that have modest optical 
depths, while the Lyman a emission becomes zero 
around line center. 

One can show that 



e x e 



(4) 



with t — t p h, as more and more photons get 
trapped in the line core for times exceeding the 
dynamical time. The multiplier (3 ~ 2 — 3 in- 
corporates details of the gravitational collapse of 
gas shells (e.g. geometry, kinematics) and does not 
impact our results as long as it does not (or only 
weakly) depend on density. 

That is, the decrease in the number of escap- 
ing/cooling Lyman a photons is approximately 
proportional to the total number of photons some- 
where in the line multiplied by the average time a 
given photon spends in the medium per unit of 
free-fall time: gravitational collapse scales with 
tff oc n -1 / 2 oc r and the number of already col- 
lapsing shells a photon would have to traverse thus 
increases linearly in space and time, i.e., -de ~ 
edt/tff. It is implicitly assumed here that the 
scattering-broadened line width, typically larger 
than 200 km/s (Dijkstra et al. 2006), exceeds any 
systematic velocity shifts, which is a good approx- 
imation for large optical depths. Hence, in Equa- 
tion (2a) the factor e is now competitive with the 
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temperature dependence of Lyman a cooling be- 
cause it picks up an exponential function of den- 
sity. 

Typically, one has tff <~ 1.6 x I0 15 /n^ 2 s and 
t ph ~ 5.0 x 10 li /n 1 /9 (M h /10 9 M Q ) 1 / 3 s for the 
adopted halo characteristics. Note in this expres- 
sion, the weak and negative dependence of t p h on 
density. This is a consequence of the random walk 
in both coordinate and frequency space that is per- 
formed by the Lyman a photon, yielding a weak 
Tq^ 3 ~ n 2 / 9 dependence for t p h, while the size of 

— 1/3 

the halo scales as n for a fixed mass M^. 

The expression for the local thermal balance, 
Equation (1) above, formally does not change, al- 
though all quantities pick up a time dependence, 
as long as local thermal balance holds. This is still 
true for t ph > tff and T ~ 10 4 K, given that the 
time needed to thermalize through collisons scales 
as l/n and the free-fall (heating) time as 1/n 1 / 2 . 
Similar considerations apply to the ionization bal- 
ance of hydrogen, but the presence of shocks would 
require a more careful treatment. The velocity 
gradients that exist maximally have a magnitude 
of 

Sv ~ r/tf f (r) ~ 10 2 km/s, (5) 

smaller than the scattering broadened line width, 
and are independent of r if an isothermal density 
distribution pertains. 

One can determine 7 straightforwardly for t ~ 
tff and a fixed mass M^. One finds that 

I + J-Bn 7 / ls 
7 ;o 5 (Cn 1 /2) + Bn 7/is' > 

where it should be noted that t p h/tff oc n 7 / 18 and 
that B w 0.5 - 0.1 cm 7 / 6 for M h = 10 9 - 10 7 M© 
(so Cn 1 / 2 « Bn 7 / 18 ). 

Evaluation of Equation (6) yields 7— 1 ~ 0.01 — 
0.5 for hydrogen densities of 1 — 10 5 cm" 3 for z = 
20 and a 10 8 M Q halo. Note that a density of 1 
cm" 3 is achieved for our halos after a contraction 
in radius by a factor of a few, much less than the 
contraction factor A -1 ~ 20 after which a disk 
forms (Mo, Mao & White 1998). One finds, for 
the 10 9 Mq halo at z = 20 and for the appropriate 
(column) density scaling with r, e.g., r <~ r" 1 
and n <~ r" 2 , that tff > t ph and 7 > 1.1 for 
T < r d V n = 0.02L. This implies enclosed masses, 
M ~ r, of about 0.02 - 0.003714 for the adopted 



isothermal profile and halo masses. The adiabatic 
value 7 = 4/3 is achieved for r < rd yn = 0.0021/ 
and a 10 9 M Q halo at z = 20 (but see see some 
corrections to 7 in the next subsection). 

The presence of the 'C term' from Equation (3) 
does not mean that conversion of Lyman a pho- 
tons to the two-photon continuum is a significant 
sink. Rather, for the considered halos, trapping of 
Lyman a occurs already at densities for which col- 
lisional de-excitation by protons is negligible (de- 
spite the large columns Nh ~ N c or somewhat 
smaller) because the thermal electron abundance 
is very small 1 . We will come back to the conse- 
quences of the rise in temperature associated with 
7 > 1 in the next subsection. Finally, the 7/18 de- 
pendence on density renders our results relatively 
insensitive to subtleties in the Lyman a radiative 
transfer. 

3.3. Dynamic Case: Two-Quantum and 
Hell Corrections 

We have assumed that the gas remains close to, 
but not exactly at, T = 10 4 K as far as its den- 
sity profile is concerned. This is reasonable given 
the sharpness of the Lyman a cooling function. A 
value 7 > 1 implies of course that the temperature 
rises with increasing density, but Lyman a cooling 
will dominate the local thermal balance for tem- 
peratures T < 5 x 10 4 K. Of course, as the tem- 
perature rises, so do the electron and proton abun- 
dance and this favors the two-photon continuum 
by decreasing N c . From Equation (6), thermal 
ionization balance and the results for N c (y), one 
finds that 7 weakens towards unity above ~ 2 x 10 4 
K when the electron abundance exceeds 0.1. How- 
ever, this temperature is reached when the density 
is 10 3 - 5 cm" 3 and two-quantum decay is quickly 
shutdown during the collapse as a density of 10 5 
cm" 3 is exceeded. 

In fact, at temperatures above 5 x 10 4 K Hell 
line cooling dominates, and the latter also suffers 
from photon trapping when hydrogen columns ex- 
ceed 10 24 cm" 2 . That is, the Hell Lyman a line at 
304A is similarly opaque (Neufeld 1990), barring 
the appropriate changes in Einstein A coefficient, 
elemental abundance etc., as its HI counterpart. 

In fact, at densities below 10 3 cm" 3 one has Nh << N c , 
but we retain the intuitive form of Equation (6) because 
the logarithm renders any error insignificant anyway. 
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One has that T He = 5.2 x 1CT 14 7V(iJe+)T 4 1/2 . 
Hence, photon trapping will continue, for the mas- 
sive halos that we consider, into the Hell regime at 
large columns. Also, the much larger Hell Lyman 
a Einstein A coefficient of ~ 10 10 s _1 boosts the 
required value of y for a given N c by two orders of 
magnitude. The Hell two-photon channel is shut 
down since densities exceed 10 5 5 cm" 3 around the 
Hell cooling peak (A 2s i s ~ 8.2Z 6 s _1 ). In any 
event, most Hell two-quantum decay photons are 
absorbed by the (HI Lyman a trapping) neutral 
hydrogen that surrounds the halo core. 

As a result of all this 7 remains well above unity 
and the system evolves adiabatically for densities 
above ~ 10 5 cm~ 3 . Equation (6) provides a good 
fit to the atomic physics of HI and Hell between 
n = 1 and n — 10 7 cm~ 3 if corrected for HI two- 
photon decay and the change in line opctical depth 
as HI cooling is superseded by Hell cooling. One 
finds that 



7-1 w-- 



^B'(n)n 7 / ls 



Zo 3 (Cn 1 /2) +B '(n)n 7 /i8' (7) 

where B' = for n = 10 3 (M, l /10 9 M Q )- 1 /3 _ 10 5 
cm~ 3 following Equation (3), and where B' w 
0.36B if n > n c = 10 5 5 cm" 3 and B' = B other- 
wise. Note here that t p h > tff always holds and 
that B'n 7 / 18 >> 1 for densities larger than n c , 
and for all halos, i.e., the HI to Hell switch has 
a modest impact because the optical depth enters 
into tph with a 1/3 power. 

4. Discussion and Future Work 

The Jeans mass for a 0.1 cm~ 3 halo is about 
Mj ~ 3 x 10 7 M Q at 10 4 K. Over the Lyman 
a cooling regime, Mj decreases only by a fac- 
tor of about 8. That is, a 0.02 - 0.003M h core 
will likely not experience significant fragmenta- 
tion during gravitational collapse up to densities 
of - n(r c ) « 10 5 cm~ 3 at z = 10 - 20, af- 
ter which fragmentation is halted adiabatically. 
That is, the system cannot cool above a few times 
10 5 K either because photons produced by, e.g., 
bremsstrahlung cannot escape since the bulk of 
these cooling photons are at energies above a few 
times 10 15 Hz and are re-proprocessed into Lyman 
alpha and trapped in the surrounding neutral ex- 
terior of the collapsing cloud. Also, the rise in 7 is 
moderate enough to justify our use of an isother- 
mal density profile. 



Although our results are order of magnitude 
estimates, they connect quite well with the de- 
tailed numerical simulations of gravitational col- 
lapse by Jappscn et al. (2005) and Klcsscn, Spaans 
& Jappsen (2005). The former authors find that a 
switch to a 7 > 1 region in density space for a col- 
lapsing gas sets a characteristic mass scale for frag- 
mentation through the Jeans mass at the ambient 
density and temperature. Hence, a value 7 > 1 
appears to be a robust indicator of the lack of 
fragmentation. As such, the picture that emerges 
from detailed hydrodynamical simulations and the 
shape of the EOS is at least consistent. 

Finally, the frequency shift that a Lyman a 
photon experiences before escape scales approx- 
imately as fshift ~ 7 11 / 6 (Dijkstra et al. 2006). 
Hence, fluctuations in temperature do not strongly 
influence our results in this respect either. 

Thus, allowing for some expected inefficiency, 
we infer that of order 0.1% of the baryon mass 
forms a pregalactic black hole of mass Mbh ~ 
10 4 - 10 6 M Q . Note here that Bromm & Loeb 
(2003) find a similar inhibition of fragmentation 
from detailed hydrodynamic simulations that as- 
sume a rougly isothermal (7 ~ 1) collapse. They 
do not include the trapping effects discussed here. 
So unless there are numerical resolution effects 
that play a role, a value of 7 = 1 may already 
be sufficient to halt fragmentation. 

In any case, these so-called intermediate mass 
black holes (IMBH) are plausible seeds for gener- 
ating the supermassive, ~ .001M spheroid, black 
holes found in galaxy cores today (c.f. Hiring 
& Rix 2004) by gas accretion (Islam, Taylor & 
Silk 2003). The inferred presence of pregalactic 
IMBH and their associated accretion luminosity 
has been a source of intense speculation with re- 
gard to a mechanism for the reionization of the 
universe (e.g., Madau et al. 2004; Rlcotti & Os- 
triker 2004; Venkatesan, Giroux & Shull 2001). 
Our results place these speculations on a sounder 
footing. Moreover, isolated IMBH should exist in 
galactic halos at a similar mass fraction accord- 
ing to simple models for generating the Magorrian 
correlation between central black hole mass and 
spheroid velocity dispersion (Islam, Taylor & Silk 
2004; Volonteri, Haardt & Madau 2003) and pos- 
sibly be detectable as gamma-ray sources (Zhao & 
Silk 2005). 

Still, there are a number of other issues that 
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should be addressed in the future. 

la) The stiffening of the polytropic equation 
of state found in this work depends crucially on 
the absence of any H2 molecules. For tempera- 
tures above 3,000 K this seems plausible because 
collisional dissociation and charge exchange with 
H + limit the abundance of H 2 , while Lyman a 
trapping keeps the temperature above 10 4 K. Still, 
H 2 may also form directly in massive halos with 
virial temperatures above 10 4 K (Oh & Haiman 
2002) and reach a universal abundance of ~ 10~ 3 . 
H 2 formation in these cases is a consequence of 
a freeze-out of the H 2 abundance, in the pres- 
ence of a large free electron fraction, as gas cools 
from above 10 4 K on a time scale that is shorter 
than the H 2 dissociation time. However, in the 
dense and massive halos that we consider Lyman 
a trapping, already at modest densities of 1 cm~ 3 , 
causes the cooling time to increase exponentially, 
from a level of a few xlO 6 yr at 8,000 K, and 
to remain larger than the H 2 dissociation time. 
That is, the gas lingers at 10 4 K, unable to reach 
8,000 K or less, and stays at those temperatures 
because H 2 formation is suppressed (H 2 is easily 
destroyed by H + ) at the ambient temperatures 
(Oh & Haiman, reaction 17 on their page 15). 
As a consequence the H 2 abundance remains at 
a very low level around 10 4 K, see Figure 4 of Oh 
& Haiman (2002), and does not contribute to the 
cooling. Furthermore, Figure 7 of Bromm & Loeb 
(2003) shows that even a halo with T vir ~ 10 4 K 
(baryonic mass of ~ 10 7 M Q ) first heats a large 
part of the cold (30-100 K) infalling gas to tem- 
peratures of 3, 000 — 10, 000 K for densities of ~ 1 
cm~ 3 , when the formation of, and cooling by, H 2 is 
incorporated. This is the relevant, minimum tem- 
perature range because there is a trough between 
the Lyman a and H 2 cooling curves here (Oh & 
Haiman 2002). In this, it is important to realize 
that our halos are quite massive, baryonic masses 
of 10 7 - 10 9 M , and dense (z > 10). These values 
favor Lyman a trapping. 

lb) Still, the presence of a UV background from 
popIII stars, that suppresses the abundance of H 2 
molecules (Bromm & Loeb 2003), would certainly 
be welcome. The critical density for collisional H 2 
dissociation to dominate is about 300 cm~ 3 , if a 
UV background as in, e.g., Bromm & Loeb (2003) 
is present with which H 2 collisional dissociation 
has to compete. We do reach this regime early in 



the collapse so that self-shielding effects would not 
limit the benefits of H 2 photo-dissociation much 
(see Bromm & Loeb 2003). Of course, in the ab- 
sence of a background radiation field, all the time 
scales (formation, dissociation, etc.) in the sys- 
tem scale as 1/density and thus their ratios are 
independent of density and the discussion of la 
applies. 

lc) The formation of the black hole will intro- 
duce a quasar whose power-law spectral energy 
distribution can boost the formation of H 2 through 
the H~ route (Haiman, Abel & Rees 2000). Hence, 
for rcdshifts below ~ 300, where H~ is no longer 
destroyed by the CMB, black hole formation as 
discussed here will facilitate the formation of H 2 
and impact the EOS of the gas surrounding the 
black hole (Scalo & Biswas 2002). The mode that 
we describe here would then be inhibited unless 
an H 2 dissociating UV background as in Bromm 
& Loeb (2003) or Oh & Haiman (2002) is present. 
Still, the large columns that we consider would 
shield at least part of the gas from X-ray feedback 
(see 2b below). 

2a) Trace amounts of dust as little as 10~ 4 of 
Solar arc sufficient to absorb all Lyman a photons 
in a homogeneous halo for the columns considered 
in this work (Neufeld 1990). Hence, the stiffening 
of the EOS that we have found disappears once 
the first metals have been produced because dust 
emission is optically thin. Also, metals are effi- 
cient coolants and, if present, would take over the 
cooling for radial contraction factors larger than 
60 (Haiman & Spaans 1999). Inhomogeneity sup- 
presses dust absorption, but facilitates the escape 
of Lyman a photons by boosting e (Haiman & 
Spaans 1999). In any case, the formation of these 
massive black holes is stopped once the ambient 
metallicity increases due to star formation. Hence, 
the fraction of massive primordial galaxies that 
harbor these black holes is dictated by the frac- 
tion of metal-free gas at z = 10 — 20. Given that 
popIII star formation is co-eval with this epoch, 
the metal-free gas fraction is uncertain and clumpy 
(Scannapieco, Schneider & Ferrara 2003). Hence, 
the overall contribution of this mode of black hole 
formation is somewhat undetermined and likely to 
lie anywhere between 5 and 50 percent, depend- 
ing on the proximity of other galaxies. If efficient, 
this mode may violate the three-year WMAP con- 
straint on the electron scattering optical depth of 
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r e « 0.09 (Spergel et al. 2006), because of the large 
X-ray output expected for these massive and late 
black holes. 

2b) Fortunately, following Ricotti et al. (2005) 
it is possible to constrain the contribution to r e . 
The latter authors find that r e scales as l/log(Nn) 
and levels off to t c ~ O.f for columns in excess of 
fO 22 cm~ 2 . Given that all our black holes form as 
the end product of a central collapse from a mas- 
sive halo, the surrounding gas has column densi- 
ties between ~ 10 22 cm~ 2 (from the initial halo 
masses and redshifts) and ~ 10 26 cm -2 (from the 
7 = 4/3 adiabatic points). We have used the mod- 
els of Meijerink & Spaans (2005) to confirm that 
these columns are sufficient to re-process X-rays 
between 1 and 30 keV for metallicities in the sur- 
rounding gas between and 10 -2 . Hence, a value 
of r e ~ 0.1 is appropriate for our black holes, even 
if they would be the dominant mode of black hole 
formation. 

3) We assume zero angular momentum for the 
gas, but the radiative transfer is not sensitive to 
the associated velocity field. Of course, our ar- 
guments are purely thermodynamic in nature and 
do not solve the angular momentum problem if the 
initial cloud is rotating. 

4) We would expect that dwarf spheroidals 
should have central black holes in the range of 
10 4 -10 6 M Q , whereas irregulars, and in particular 
late-forming dwarfs, should not have such central 
IMBH. In order to substantiate this, more detailed 
hydrodynamical simulations that include dynami- 
cal photon trapping, or its EOS parameterization, 
should be performed. 
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